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Abstract 



Actin networks, acting as an engine pushing against an external load, are 
fundamentally important to cell motility. A measure of the effectiveness of 
an engine is the velocity the engine is able to produce at a given force, the 
force- velocity curve. One type of force- velocity curve, consisting of a concave 
region where velocity is insensitive to increasing force followed by a decrease 
in velocity, is indicative of an adaptive response. In contrast, an engine 
whose velocity rapidly decays as a convex curve in response to increasing 
force would indicate a lack of adaptive response. Even taken outside of a 
cellular context, branching actin networks have been observed to exhibit 
both concave and convex force- velocity curves. The exact mechanism that 
can explain both force-velocity curves is not yet known. We carried out an 
agent-based stochastic simulation to explore such a mechanism. Our results 
suggest that upon loading, branching actin networks are capable of remod- 
eling by increasing the number filaments growing against the load. Our 
model provides a mechanism that can account for both convex and concave 
force-velocity relationships observed in branching actin networks. Finally, 
our model gives a potential explanation to the experimentally observed force 
history dependence for actin network velocity. 

Key words: Branching actin networks; actin-based motility; force- velocity 
relationship; stochastic simulations 
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Introduction 

Branching actin networks are the principle engine that drives cell motility 
ranging from cell migration (1, 2) to endomembrane trafficking (3). In the 
lamellipodium of migrating cells, actin filaments grow from their barbed- 
ends, pushing against the plasma membrane in the direction of cell move- 
ment. New filaments branch off of existing filaments through the actin 
related protein (Arp2/3) complex, activated by WASP at the membrane. 
Filaments branch at a characteristic angle of ~ 70° . Capping proteins limit 
the growth of filaments by binding to the barbed-end of the filament. At the 
back of the network, actin filaments depolymerize and are severed, providing 
a fresh actin monomer supply to the front (1, 2). Understanding the basic 
process by which an actin network is able to exert force against a load is a 
fundamental step to understanding a number of cellular processes (4). 

Both in vitro and in vivo experiments have been performed to probe 
the force- velocity relationship of growing actin networks (5-10) . Migrating 
cells show an adaptive response exhibiting a concave force-velocity relation- 
ship (6, 7). However, the concave force- velocity relationship is preceded by 
a large reduction in velocity to small forces (6). The mechanism control- 
ling the concave force-velocity relationship and the initial response to small 
forces in cells is complicated by other cellular components such as focal ad- 
hesions. To study the exact mechanism that determines the force-velocity 
relationship of a branching actin network, in vitro experiments with more 
controllable conditions have been performed. One study measured the ve- 
locity of an actin network growing against a constant load force set by an 
atomic force microscope, and the resulting force-velocity relationship was 
convex (10). In a different in vitro experiment, the actin network grew 
against the flexible cantilever. The load force thus progressively increased 
as actin polymerization deflected the cantilever, and the network showed 
a concave force- velocity relationship (9). That experiment also showed a 
hysteresis effect where the velocity of the network was dependent upon the 
past forces applied to the network. 

Experiments done in vitro have demonstrated both convex and concave 
force-velocity relationships in branching actin networks. This suggests that 
actin networks can respond to external forces in both an adaptive and a 
non- adaptive manner outside of cellular context. Even within the simplified 
in vitro setting, it is still unclear how the individual factors that govern 
branching actin network dynamics generate both the convex and the concave 
curves. It has previously been proposed that the actin network remodels 
itself in response to force (9), but the nature of such remodeling is largely 
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unknown. 

Evidence suggests that actin filaments utilize thermodynamic free en- 
ergy to add additional monomers to exert force towards the leading edge of 
the network (11). The proposed model for that behavior has been termed 
the Brownian ratchet (12-14). The Brownian ratchet mechanism takes ad- 
vantage of the asymmetry in the on and off rates of actin monomer binding 
to an existing filament. Small gaps arise between the actin filaments and 
the leading edge due to thermodynamic fluctuations. Monomers are able 
to bind during such fluctuations and push the leading edge forward. The 
predicted force-velocity curve for a single filament is a convex negative ex- 
ponential function. When many filaments grow against the same load, they 
are able to share the load. A simple force-sharing mechanism predicted that 
the force- velocity curve is nonetheless convex (15), even though its slope 
is shallower than that for single filament. A model that tries to explain 
the force-insensitive region of the concave force- velocity curve (16) is largely 
incapable of accounting for the convex curve. 

Our theoretical model was built to study both types of force-velocity re- 
lationships for branching actin networks. We used an agent-based stochastic 
simulation method inspired by Wcichsel et al. (17) and Schaus and Borisy 
(18). Our results show that the balance between growth, branching and 
capping events controls the ability of a branching actin network to rein- 
force itself against a load. The model can explain both convex and concave 
force- velocity relationships. 

Model Description 

To discern the physical mechanism governing the force- velocity relationships 
of branching actin networks, we aimed to construct the simplest model able 
to reproduce the observed effects without compromising the physical essence 
of actin networks. The model therefore only includes the four essential 
processes of branching actin networks in in vitro conditions. Below, we 
describe the qualitative features of our model. 

1. Filaments grow by adding new monomers to the barbed end of the 
filaments. When the filament is not in contact with the load, it does not 
feel the load and, hence, it could grow at its free rate Vq. When a filament 
is in contact with the load surface, the rate of growth follows the Brownian 
ratchet mechanism, where the growth velocity of the filament is reduced by 
a Boltzmann factor: 

(1) 



B{F) = exp 



FS 
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where F is the force felt by the filament, S is the length of an individual actin 
monomer, /c^ is Boltzmanns constant, and T is the absolute temperature. 

2. New filaments are created by branching off existing filaments. Arp2/3 
binds to existing filaments and creates a site for a new filament to grow and 
generates the characteristic 70° angle in between the original and newly 
branched filaments (19). Branching was modeled as a zero-order reaction, 
independent of the number of actin filaments, which is consistent with the 
experiments that suggest WASP/Scar-mediated Arp2/3 activation is the 
limiting factor of network growth (20). And the branching angle followed a 
Gaussian distribution that centers at 70°, with a standard deviation of 5° 
(19). 

3. Capping proteins can bind to the tips of actin filaments, preventing 
them from further elongation. In in vitro conditions, the lifetime of barbed 
end-bound capping proteins is ~30 min (21). This feature is modeled by 
filament growth stopped once capped. We modeled filament capping as a 
first-order reaction: the capping rate was proportional to the number of free 
barbed ends, which is agreement with in vitro experiment measurements 
(22, 23). 

4. A significant factor in the efficiency of a growing actin network against 
a load is the ability of the network to share the load across multiple filaments 
at the leading edge, which has been a recent topic of study (15). We imple- 
mented a similar load-sharing scheme to Schaus and Borisy (15) among the 
filaments in contact with the load surface, i.e., the sum of the load force felt 
by each contacting filament is in balance with the total load. Note that the 
load forces felt by individual contacting filaments across the leading edge 
are not equal, because the closer the filament orientation is to the normal 
of the load surface, the larger share of the load force this filament will feel. 

At each time step, calculations were done in the following order. First, 
the location of the leading edge was determined by the location of the fore- 
most filament. Then, filaments in contact with the leading edge (see methods 
for definitions) were located. Filaments were next capped and branched, the 
rates of which were calculated based on the location of the filaments and 
Poisson statistics. Finally, the positions of all the filament tips were up- 
dated by the growth rate for each filament in accordance to force-sharing 
mechanism. 

Model parameters were estimated based on experiment evidence where 
possible, and they are listed in Table 1 with references. 
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Methods 

Actin Model 

Each filament was modeled as a point in a two-dimensional plane represent- 
ing the filament barbed end. The plane was bounded by a hard leading 
edge in the principle direction of motion (X) and periodic boundaries in 
the perpendicular (Y) direction. Each filament had three properties: an 
X-coordinate, a Y-coordinatc and an angle of growth 9 which was relative 
to the X-direction. As the filament grew, the X and Y coordinates would 
change in time, e.g. xj = Vj cos(^j) for filament j growing at an angle 6j, 
while the angle (9) for each filament did not change. The filaments were 
limited to having angles in —90° < 9 < 90° as filaments growing against 
the principle direction of growth (X) are not seen in experiment (24) and 
would quickly grow too far away from the leading edge to contribute to the 
network velocity. The default parameters for the model are listed in Table 
1. 

The filaments grew against the force exerted by the leading edge (at 
point X=0). A filament was treated as being in contact with the boundary 
if the tip of the filament was within one subunit length (5) from the leading 
edge. Identifying the filaments that were in contact with the leading edge 
was performed at the beginning of each time step. 

Within the simulations, new filaments were generated by branching from 
existing filaments. The filaments branched at a constant rate, calculated us- 
ing zeroth-ordcr Poisson statistics. Having a constant rate of branching new 
filaments essentially assumes that the rate limiting factor is the concentra- 
tion of Arp 2/3 (16, 17, 20, 25). The initial tip of a new branch was placed 
at a point along the initial filament randomly selected from a uniform dis- 
tribution extending a distance oi 5 5 backwards from the tip of the initial 
filament. The difference between the angle of the branching filament (^5) 
and the angle of the initial filament {9i) was drawn from a Guassian distri- 
bution (A/") with mean 70° and variance cr^ = 25°, i.e. {9h - 9i) ~ Ar(70, 25), 
and the filaments branched in both directions (^5 > 9i or ^5 < 9i) with equal 
probability (19). Gaussian random numbers were drawn using the gasdev 
algorithm in (26). The filaments were capped as a first-order reaction, with 
first-order Poisson statistics, meaning the capping rate was dependent upon 
the number of active filaments not in contact with the leading edge. Once 
a filament was capped, it was no longer able to grow, branch or exert force 
upon the leading edge (21). The first-order statistics were calculated based 
on the total number of active, uncapped filaments not in contact with the 
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leading edge. At each time step, the number of filaments added and capped 
during the time were calculated using a uniform (0,1) random number using 
the rani algorithm in Numerical Recipes (26). Filaments in contact with 
the leading edge were neither branched nor capped as in previous models 



Filaments grew at a deterministic, constant rate in the direction defined 
by 9, based on the assumption that the rate of actin subunits adding to 
the filament is much faster than the rate of the network reorganization. We 
were principally interested in the geometric structure of the network and 
considered only average growth dynamics. Similar to in vitro conditions, 
the actin monomer concentration was modeled as saturated leading to a 
constant growth rate of 100 6 /s. The velocity of each individual filament in 
contact with the leading edge was reduced by the Boltzmann factor derived 
from the Brownian-ratchet mechanism (12-14). The force was shared across 
all filaments using a modified version of the optimal force sharing from 
Schaus and Borisy (15). Each time a filament would add an additional actin 
subunit, the energy penalty would be proportional to FtotScos{9), where F 
is the total force applied by the leading edge. Normalizing the individual 
forces so that the total force adds up to Ffot gives the final velocity expression 
for a filament i in contact with the leading edge: 



where the sum j is over all of the filaments in contact with the leading edge 
and Ftot is the total force applied to the system. 

At each time step, calculations were done in the following order. First, 
the location of the leading edge was calculated, which was defined as the 
largest X coordinate of the active filaments. The filaments which were then 
within S of the leading edge are marked as in contact with the leading edge 
and were no longer able to be capped or branch new filaments. Next, the 
number of branching and capping events between times t and t + dt was cal- 
culated by comparing a uniform random number in (0,1) and the cumulative 
distribution function for the associated Poisson distribution. Individual fila- 
ments were chosen randomly from the population of filaments not in contact 
with the leading edge to be capped. Following that, filaments were again 
randomly chosen to serve as the branching point for new filaments. Next, 
the normalizing constant was calculated by summing over the filaments in 
contact with the boundary to get the average force exerted by each filament. 
Finally, the positions of the filaments are advanced by Vidt where Vi is the 



(15, 18, 27). 
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velocity of filament i. 

All simulations started with independent, identically distributed initial 
conditions. The simulations started with 200 filaments uniformly distributed 
in the y direction and uniformly distributed within a box of length 20 6 in 
the X direction with the constraint that half of the filaments were in the 
first 10 5 and half in the second 10 5. The initial filament orientations were 
randomly drawn from a uniform distribution on (—90°, 90°). The time step, 
dt, used in each simulation was 10~2s. 

Short Time-Scale 

One thousand simulations were run with independent, identically distributed 
initial conditions as above. The simulations were averaged at each time step 
to generate Figure 1(A). Figure 1(B) was generated by taking the minimum 
velocity after force was applied. Parameter values were chosen to be At = 1 
and A = 200 to emphasize that we were only able to see convex force velocity 
curves on short time-scales. 

Equilibrium Simulations 

Each data point from Figure 2 was the result of averaging 10 equilibrium 
simulations with constant force starting from the initial conditions described 
above. Simulations were run for a total of 10,000 s with a dt of 10~^ s and 
data was sampled every 0.2 s. Reported data was sampled from the second 
half of the simulation to allow the system to minimize the influence of initial 
conditions. The equilibration can be seen in supplemental Figures Sl-2. 

Hysteresis Simulations 

The hysteresis simulations were performed with the same initial conditions. 
Simulations were run for a total of 7,500 s. The first 2,500 s were run at 
an initial (low) force /q to equilibrate the system to the velocity observed 
in Figure 2. At time 2,500 s, the force was increased to /i. The system 
was allowed to equilibrate again to a lower velocity. Finally, at time 5,000 
s the force was reduced back to /o where the velocity rapidly rose before 
converging back to equilibrium. 

Filaments were made to 'stick' to the leading edge by reducing the 
amount of force felt by a filament close to 6 away from the leading edge. 
The new expression force felt by each filament was Fj = F'^q{x) where F* 
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is the force used in (2) and q{x) is a cubic polynomial such that: 

{1 if < a; < 7 

-2[^y +3[^y if7<x<i (3) 
if X > 1 

where 7 = 0.9 for our simulations. 

Curve Pcirameter 

We characterized the force- velocity curves by their /i/2) the force at which 
the relative velocity was reduced to ^. Twenty constant force simulations 
were performed for each pair of branching and capping rates ranging from 
850 pN to 17 nN. fi/2 was estimated by linearly interpolating between the 
two successive forces where the velocities surrounded ^ . When the velocity 
at 17 nN was greater than one half, we used the largest force calculated (~17 
nN) for the purposes of figure 5. Likewise, for the purely convex curves with 
/1/2 < 4 nN the value 4 nN was used in the graph for clarity. 

Results 

Our simulations were focused on how the collective properties of a branching 
actin network influence the ability of the network to grow against a flat 
surface applying a load force. 

Fast force-velocity response is always convex 

The first set of simulations we ran tested the temporal response of the ve- 
locity of branching actin network against a fixed load force. Figure 1(A) 
shows that upon loading force, the velocity drops almost instantly, which 
then recovers in a longer timescale (~minutes) reaching a value lower than 
before the application of a load. Taking the velocity at the bottom of the 
initial response to force response gives a force-velocity curve for the short 
time-response of the network. Our model suggests that the force-velocity 
curve is always convex (Figure 1(B)). 

Force- velocity relationship at long-time response can be either 
convex or concave 

Running the simulations for an extended amount of time allows us to study 
the equilibrium force velocity relationship. Figure 2(A) shows that we were 
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able to reproduce both convex and concave force-velocity curves. The only 
difference between the sets of simulations is the absolute value of the capping 
and branching coefficients. Their ratio, and therefore the average number 
of filaments, was fixed. We hypothesize that the network is able to reinforce 
itself by bringing more filaments to the leading edge. 

Branching actin networks remodel by increasing the number 
of filaments contacting the leading edge 

An actin network can be visualized as a population of filaments contacting 
the boundary and another population of filaments trailing the leading edge 
in reserve. The network reinforces the filaments at the leading edge when 
trailing filaments grow to reach the leading edge. That remodeling response 
simply depends on the rate at which trailing filaments are able to catch 
up to the leading edge. Figure 3 is an explanatory diagram showing a 
hypothetical branching pattern. The first filament is in contact with the 
boundary, the second filament branches off the first one, is further back, 
and serves as a source for new filaments. We term such filaments reserve 
filaments (Figure 3(A)). Since these reserved filaments are not in contact 
with the load surface, they dont feel the load and grow at their free rate 
faster than the leading edge. Consequently, some of the reserved filaments 
catch up to the leading edge. When the capping rate is high, these reserved 
filaments will get capped and stop growing before reaching the load surface 
(Figure 3(B)). Conversely, when capping rate is low, these reserve filaments 
grow to be in contact with the load surface, thereby reinforcing the leading 
edge against the load (Figure 3(C)). 

The hypothesis that the network reinforces itself by filaments growing to 
the leading edge would suggest that a network composed of shorter filaments, 
being less likely to grow to the leading edge, would stall at lower forces. 
Increasing the capping rate causes the filaments to grow for shorter periods 
of time leading to shorter filaments, and indeed. Figure 2(B) shows that 
by changing the capping rate, the force-velocity curve continuously changes 
from concave to convex shapes. Fewer filaments reinforce the leading edge 
leading to lower and lower stall forces. The number of filaments in contact 
with the leading edge attests the level of reinforcement. Figure 2(C) shows 
that more filaments are recruited to the leading edge when the capping rate 
is low for the same level of force. 

Figure 2(C) shows that the number of filaments in contact with the lead- 
ing edge increases with increasing force for all cases. However, the cases with 
lower capping rates generate a larger increase before peaking and stalling. 
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This is because the branching rate is zeroth order, where the rate is inde- 
pendent of the number of filaments; whereas the capping rate is first order, 
whose rate increases with the number of available filaments. Before reach- 
ing the peak, the network is, in some sense, branching dominated. On the 
other hand, increases in force increase the number of filaments in contact 
with the leading edge, so does the capping rate. When the effective capping 
rate becomes faster than that of branching, the contacting filaments will get 
capped and, hence, their number will drop off. It is also of note that the 
specific filaments in contact with the leading edge are constantly changing. 
When the capping effect dominates over branching, the leading edge velocity 
is no longer sufficient to provide the turnover to sustain the higher number 
of filaments. Approximately at the peak of Figure 2(C), the network begins 
to stall in the force- velocity curves in Figure 2B. 

Branching actin network remodeling could account for the 

observed hysteresis effect 

The accTinuilation of filaments at the boundary is also able to explain some 
of the hysteresis observed in actin networks. As seen in Figure 2, the number 
of filaments in contact with the boundary increases with increasing force. 
Subjecting the network to a large force and subsequently releasing that force 
should leave excess filaments in contact with the leading edge. 

Figure 4(A) shows the velocity of a simulation where the network pushed 
against low force for the first third of the run, followed by high force in the 
middle third, and finally the original low force. That is comparable to the 
experiment in Parekh et al. (9). The velocity initially shoots back up 
in response to the reduced force, but it rapidly decays back to the initial 
equilibrium force similar to (16). The reduction in the load force speeds 
up the growth rate of all the contacting filaments at the leading edge. Due 
to the angle dependence of the load sharing for individual filaments, the 
speed-up of growth rates is heterogeneous across the contacting filaments. 
As a result, some filaments grow faster while staying in con- tact with the 
load, and slower growing filaments slide of the leading edge and are capped. 
Ultimately, the number of contacting filaments relaxes back to the velocity 
corresponding to the original force and completely loses its memory of the 
previous loading force. Figure 4(B) shows that sustained hysteresis can 
be realized in the model, if we incorporate a factor that causes the actin 
filaments to stick to the leading edge. While we do not know the exact nature 
of such an interaction between the filament tips and the load surface, actin 
tethering has been theoretically proposed (14) and has some experimental 
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evidence (28-32). 

The balance between capping and branching events dictates 
actin network remodehng 

Our results suggest that the actin network remodels itself by changing the 
number of filaments in contact with the leading edge. That remodeling in 
turn determines the shape of the force velocity relationship. In particular, it 
determines the length of the concave portion of the curve. The balance be- 
tween the branching and capping rates controls the nature of the remodeling. 
To obtain a systematic, systematic understanding of the rate dependence of 
the actin force-velocity relationship, we carried out a phase diagram study 
for the capping and branching rates. We used the force at which the velocity 
drops to 50% of the small load force to characterize the shape of the force 
velocity curve. Figure 5 shows the estimated fi/2 values at each parameter 
value that constitutes our principle prediction: faster capping rates lead to 
less concave force velocity curves, and faster branching rates lead to more 
concave curves. 

Discussion 

We have proposed a simple mechanism where branching actin networks re- 
model against a load force. The model shows that the initial response of 
branching actin networks to loading always gives a convex force- velocity re- 
lationship (Figure 1). On longer time scales, smaller capping rates and larger 
branching rates generate more concave force-velocity relationship (Figure 2 
and 5). 

A number of recent theoretical studies have focused on how to explain 
both the convex and concave force- velocity relationships for branching actin 
networks (17. 33). Likewise, the nature of hysteresis effects observed in 
experiment (9) remains a subject of inquiry. Multiple attempts have been 
made to explain the stall force of individual actin filaments, but the stall 
force of a network of cooperating actin filaments is poorly understood (24) . 
Our simulations yielded a stall force of approximately 2-3 pN/filament (see 
supplement), which is in close agreement with the reported value of 1.7ib0.8 
pN/filament (7). Our reported stall force per filament provides evidence 
that actin networks use close to optimal force sharing. 

The importance of the number of growing filaments at the boundary 
determining how the network responds to load force has previously been 
suggested for branching actin networks (9, 16, 27), and for bundled filaments 
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such as actin (34) and microtubules (35, 36). In the context of our model, 
bundled filaments (actin or microtubule) qualitatively correspond to the 
case of a very low branching rate or a high capping rate. According to our 
calculation, although the number of contacting filament in this case increases 
with the load force, the reinforcement is limited (the lower two curves in 
Figure 2C). Consequently, the resulting force- velocity curve is always convex, 
consistent with the previous findings (34-36). 

Our model therefore suggests that branching event is essential in yielding 
a concave force-velocity relationship of actin network growth against load 
(Figure 5). It should be noted that a similar conclusion was also reached 
by Carlsson's model (16). Carlsson (16) suggested that actin networks with 
autocatalytic branching would continually increase the number of filaments 
at the boundary leading to force-independent velocities. Although branching 
events are autocatalytic as demonstrated in experiments (37), Carlsson's 
model (16) by itself can only account for the force-insensitive region of force- 
velocity curves, i.e., an additional negative feedback would be necessary to 
limit the increase in density predicted by the model to produce the inevitably 
reduction in velocity at high forces. Moreover, the Carlsson model (16) 
predicts a transient hysteresis effect that does not correctly reproduce the 
sustained hysteresis observed in experiment (9). 

In contrast to Carlsson's model, our model has a built-in negative feed- 
back mechanism. The capping rate is a first-order reaction of the number 
of free filament barbed ends while the branching rate is constant, indepen- 
dent of the number of filaments. As the number filaments increases with 
the load force, so does the capping rate. That effect limits the total number 
of filaments. As such, our model explains both the concave and the con- 
vex force-velocity curves without resorting to additional mechanism (Figure 
2). Also, our model can explain the sustained hysteresis in force- velocity 
relationship (9) (Figure 4(C)). 

Our model does not preclude any other negative feedback mechanism 
limiting the density of actin filaments. It is likely that increased filament 
density would lead to excluded- volume effects at large forces (33). Any 
external negative feedback mechanism would limit the length of the force- 
insensitive region of the force- velocity curve. It is important to note that, 
due to the exponential term in Eq. (1), even a relatively small change in 
filament density would lead to a large change in velocity. A doubling of the 

r F6 1 

number offilaments, N , would lead to an approximately exp ^ - fold 

increase in the velocity. A surprising feature we observed in our simulation 
was a substantial reduction in network velocity to extremely small forces 
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(<200 pN, see supplemental figure S3). This reproduces the effect seen in 
(9), which cannot be reproduced by (33). We hypothesize that network 
velocity has to be reduced by a sufficient amount before trailing filaments 
are able to catch up to the leading edge. When the leading edge is moving 
close to the growth rate of individual filaments, trailing filaments are unable 
to reach the leading edge. When the leading edge is sufficiently slowed by 
the opposing force, trailing filaments are able to reach the leading edge and 
the leading edge velocity stabilizes. 

It is generally believed that new filament branching in actin networks 
occurs in a narrow zone near the membrane (20, 38), although there is no 
quantitative measurement on the exact location and the size of such active 
Arp2/3 complex-enriched zone. It has been suggested that the restriction 
of filament branching at the membrane may have a role in the geometric 
organization of the network (18, 38). We ran simulations to test if restricting 
the area where new filaments could branch would change the predictions of 
our simulations. We restricted the branching of new filaments to a zone 
of distance L away from, but not at the leading edge, while keeping all 
the other model parameters the same as those in Figure 2(B). Here, the 
distance L ranges from 2 to 40 actin subunit lengths, corresponding to 5.4 
- 108nm. Figure S4 shows that the spatial restriction of branching events 
(L = 5.4nm - lOSnm) does not significantly change the qualitative behavior 
of force-velocity relationship as compared to the unrestricted case (Figure 
2(B)). If we allowed the branching events strictly at the leading edge, all the 
resulting force- velocity curves become convex (Figure S5). As experiments 
show that the Arp2/3 complex itself constitutes the first subunit of the 
daughter branch (39), the branching point should be at least one subunit 
length away from the leading edge when the first actin subunit adds to the 
bound Arp2/3 complex. We therefore deem this case of branching events 
strictly at the leading edge may not reflect the reality. Nonetheless, Figure 
S4 and S5 suggest that regulation of the active zone of Arp2/3 complex 
could modulate the quantitative behavior of the force- velocity relationship 
of branching actin network. 

The majority of our simulations were performed with constant force. 
However, our model is relevant to both constant and non-constant force 
because it requires no equilibrium assumptions. Our model showed a re- 
laxation time (^minutes) before the network reached equilibrium velocity. 
Changing the force more slowly than this relaxation time would allow the 
network to continuously adapt to the increased forces and show strong hys- 
teresis effects. Increasing the force substantially faster than the relaxation 
time would not allow the network to restructure itself leading to constant 
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force type results; such phenomena are captured by the convex force-velocity 
relationship predicted by our model (see Figure 1(B)). An experimentally 
observed value for this relaxation time could be found, and experiments 
changing the force pressing against an actin network more slowly than the 
observed relaxation time could test our model prediction. 

Figure 2(B) should be qualitatively reproducible in experiment and could 
serve as an excellent test of several hypotheses for branching actin networks. 
The relevant biochemical quantities could be manipulated in the vein of 
Cameron et al. (40). Such an experiment could be performed both with 
constant force (10) and with a constantly increasing force (9). We further 
note that the number of active actin nucleation factor (such as WASP) at 
the load surface can be negatively regulated by the number of free barbed 
ends at the leading edge (41). This additional negative feedback mechanism 
could limit the effect of increasing the Arp2/3 concentration. Interestingly, 
capping proteins could increase the growth rate of branching actin networks 
by promoting more frequent filament nucleation by Arp2/3, funneling actin 
monomers to the uncapped barbed ends of actin filaments, without affecting 
the free filament elongation rate (41). These seminal in vitro studies point 
to a more intertwined interaction between branching and capping events of 
actin networks, which will be the future extension of our current model. 

Our model implements a load-sharing mechanism where the contacting 
filaments collectively share the load across the leading edge. That is, the 
addition of a new actin monomer is only opposed by a fraction of the total 
load force pressing on the network. In the context of our model, load-sharing 
mechanism is valid as long as Brownian ratchet mechanism holds up. The 
Brownian ratchet mechanism assumes that the thermal fluctviations between 
the filament tip and the load surface are significantly faster than the addition 
of new actin monomers (13, 14). Fluctuations must be large enough for a 
new monomer to fit in the gap between the tip and the load. Smaller and or 
slower fluctuations would then reduce the efficiency of the mechanism. Ex- 
periments have indeed demonstrated that reducing thermal fluctuations by 
lowering the temperature strongly hinders the efficiency of filament growth 
(11). Thus, it is the thermal fluctuations that buffer between the contacting 
filament tips and the load surface, providing the flexible interface to accom- 
modate insertions of actin monomers. The separation of time scales also 
implies that the load force felt by each contacting filament is an average 
over many fluctuations. Consequently, only the partial load force shared 
across filaments dictates the network growth rate. In the future, we intend 
to study how thermal fluctuations influence the efficiency of load-sharing 
mechanism in further detail. 
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In the model, the actin network was assumed to be a rigid structure 
so that does not buckle nor break down. In reality, there will be many 
cross-linker proteins that stiffen the actin network. In addition, capping 
events in our model limited filaments to an average length of less than 1 
/xm, significantly less than the persistence length of a single actin filament, 
~17 //m (42). Consequently, each individual actin filament can be viewed 
as a rigid rod. The actin filaments in our case are highly branched, which is 
believed to be much more rigid than its unbranched counterpart. Although 
we did not explicitly incorporate these known properties, our model used 
their effects and simply assumed that the network was a rigid structure. 
The model does not consider the cases where the filaments could undergo 
buckling or even breaking down, the topic of which will be investigated in 
the future. 

Conclusion 

The simple physical model shown here gives insight into the behavior of 
branching actin network remodeling in the presence of a load. In particu- 
lar, the network velocity dependence upon the number of filaments growing 
against the leading edge provides a simple mechanical mechanism to explain 
a number of experimental effects. The ability of actin networks to remodel 
is controlled by the balance between branching and capping rates. This 
mechanism can account for both the observed convex and concave force- 
velocity relationships. Further investigation into actin network properties, 
both physical and biochemical, that determine how many growing filaments 
a network is able to recruit to the leading edge will deepen our understanding 
of actinbased motility. 
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Table 1: Table of model parameters 



S 2.7 nm Length of an actin subunit. (43) 

ksT 4.6 pN-nm Absolute temperature. 

A 200 — 1, 600/s Branching rates. -"^ 

K 1-8/s/filament Capping rates. (21) 

Vfree 100 (5/s Default velocity. (43) 

Nf 200 Average number of free filament barbed ends. 

= A/«. 

Qbr 70° Mean branching angle. (19) 

(Tftr 5° Branching standard deviation (19) 



^Branching rate was set relative to the capping rate to determine the average 
number of free filaments. 
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Figure Legends 

Figure 1 

A) shows a characteristic time trace of the velocity response to force apphed 
at 100 seconds, and B) shows the convex force-velocity curve generated by 
the initial response of the network to force. These simulations were run with 
capping rate k = 1/ s/ filament and branching rate A = 200/s. 

Figure 2 

Varying the capping rate (k) over an order of magnitude changes the shape 
of the force-velocity curve. A) shows that we are able to qualitatively re- 
produce both the convex and the concave force velocity curves. B) shows 
the continuous deformation of the force-velocity curve for a few capping 
rates. C) shows the relationship between force and the equilibrium number 
of contacts for the same capping rates as in B. The error bars in B) and C) 
represent the standard deviation estimated from 10 simulations. 

Figure 3 

This diagram shows how the average length of the filaments could influence 
the number of filaments in contact with the leading edge. Filaments are 
represented by black lines; the red line is the leading edge; the blue circles 
are actively growing barbed ends; and the yellow circle represents a barbed 
end that has been capped. A) A hypothetical scenario involving one filament 
barbed contacting the leading edge with a second filament growing behind 
the leading edge. B) When the capping rate is high, the reserve filament 
is capped (yellow circle) before it reaches the leading edge. C) When the 
capping rate is low, sufficiently long filaments can grow and contact the 
leading edge, increasing the leading edge velocity. 

Figure 4 

These are plots of simulations testing for hysteresis in the force-velocity 
relationship. A) shows a typical simulation result with a transient hysteresis 
effect. B) shows the time dependent force used in both A and C. C) shows 
a typical simulation when filaments reaching the leading edge stuck to the 
leading edge. 
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Figure 5 

The color represents the estimated /i for a variety of branching and capping 

2 

rates. The cases with concave force-velocity curves are labeled with white 
asterisks. Decreasing the capping rate and increasing the branching rate 
serve to generate more filaments which shifts /i rightward, meaning a more 

2 

concave curve. 
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1 Equilibration 

The simulations did appear to be sampling from an equilibrium distribution. 
For the simulations where the force- velocity relationship was loosely flat, the 
velocity converged to an equilibrium quite rapidly, as can be seen in Figure SI. 
The curves shown are the average of 10 simulations with the highlighted region 
representing the estimated ± standard deviation. 
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Figure SI: A) k = 1/s/filament and F = 2.55 nN and B) k = 3/s/filament and 
F = 16.2 nN 



However, where the force- velocity curve was sharp, the velocity converged only 
at non-physiological time-scales, if at all. This can be seen in Figure S2. 



2 Stall Force 

We estimated the stall force per filament by defining the system to be stalled 
when an increase in force of 170 pN led to a decrease in velocity of less than 
1%. The stall forces per filament we observed ranged between 0.9 pN and 2.4 
pN for 1 < K < 5/s/filament with a mean of 1.3 pN. 
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Figure S2: n = 3/s/filament for both and A) F = 9.4 and B) F = 11.1 nN 



3 Velocity Reduction to Small Forces 

Our model qualitatively reproduces the large velocity reduction in response to 
small forces observed in experiment [2] . The graph below shows the equilibrium 
velocity "/^/^ee in response to a 170 pN force. We believe that this is due to the 
fact that the leading edge must be sufficiently slowed down for trailing filaments 
to catch up to the leading edge. Thus, the force-independent velocities must be 
significantly slower than w/ree- 
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Figure S3: Equilibrium velocity reduction in response to a small (^170 pN) 
force by capping rate. 
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4 Spatial Restriction of Branching 



There is experimental evidence that new filaments only branch in a small zone 
bordering the membrane. To test how well our results would hold up under that 
type of condition, we ran another set of simulations restricting where filaments 
could branch. For Figure S4, we kept all of the conditions identical to simula- 
tions in the main text except for that filaments only branched in a zone of N5 
away from the leading edge. 
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Figure S4: Force- velocity relationship for k = 3/ s/ filament where branching is 
restricted to N5. 

We also tried allowing filaments to branch directly at the leading edge. This is 
not necessarily physical as filament tips appear to nucleate new branching sites 
[1]. The results we found are not similar to current experiments as can be seen 
in Figure S5. 
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Figure S5: Force- velocity relationship for k = \/ s/ filament where branching is 
allowed at the leading edge and restricted to N6. 
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